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Dirac points lie at the heart of many fascinat- 
ing phenomena in condensed matter physics, from 
massless electrons in graphene to the emergence 
of conducting edge states in topological insulators 
[1, 2j. At a Dirac point, two energy bands inter- 
sect linearly and the particles behave as relativis- 
tic Dirac fermions. In solids, the rigid structure 
of the material sets the mass and velocity of the 
particles, as well as their interactions. A differ- 
ent, highly flexible approach is to create model 
systems using fermionic atoms trapped in the 
periodic potential of interfering laser beams, a 
method which so far has only been applied to ex- 
plore simple lattice structures [3^ Here we 
report on the creation of Dirac points with ad- 
justable properties in a tunable honeycomb opti- 
cal lattice. Using momentum-resolved interband 
transitions, we observe a minimum band gap in- 
side the Brillouin zone at the position of the Dirac 
points. We exploit the unique tunability of our 
lattice potential to adjust the efl*ective mass of 
the Dirac fermions by breaking inversion symme- 
try. Moreover, changing the lattice anisotropy 
allows us to move the position of the Dirac points 
inside the Brillouin zone. When increasing the 
anisotropy beyond a critical limit, the two Dirac 
points merge and annihilate each other — a situ- 
ation which has recently attracted considerable 
theoretical interest [5-9j, but seems extremely 
challenging to observe in solids [lOj. We map 
out this topological transition in lattice param- 
eter space and flnd excellent agreement with ah 
initio calculations. Our results not only pave the 
way to model materials where the topology of the 
band structure plays a crucial role, but also pro- 
vide an avenue to explore many-body phases re- 
sulting from the interplay of complex lattice ge- 
ometries with interactions [11, 12j. 

Ultracold Fermi gases have emerged as a versatile tool 
to simulate condensed matter phenomena O [H [13] . For 
example, the control of interactions in optical lattices has 
lead to the observation of Mott insulating phases [HI [15] , 
providing new access to the physics of strongly correlated 
materials. However, the topology of the band struc- 
ture is equally important for the properties of a solid. 
A prime example is the honeycomb lattice of graphene, 
where the presence of topological defects in momentum 



space - the Dirac points - leads to extraordinary trans- 
port properties, even in the absence of interactions [P. In 
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FIG. 1: Optical lattice with adjustable geometry, a, 

Three retro-reflected laser beams of wavelength A = 1064 nm 
create the two-dimensional lattice potential of equation ([T|). 
X and Y interfere and produce a chequerboard pattern, while 
X creates an independent standing wave. Their relative po- 
sition is controlled by the detuning 5. b, Different lattice 
potentials can be realised depending on the intensities of the 
lattice beams, as displayed above. The diagram below shows 
the accessible lattice geometries as a function of the lattice 
depths Vx and Vx- The transition between triangular (T.) 
and dimer (D.) lattices is indicated by a dotted line. When 
crossing the dashed line into the honeycomb (H.c.) regime, 
Dirac points appear. The limit ^ Vx^y corresponds to 
weakly coupled one-dimensional chains (ID c). c, The real 
space potential of the honeycomb lattice has a 2-site unit cell 
[A^B sites) and the primitive lattice vectors are perpendicu- 
lar, d, Sketch of the first and second Brillouin zones (B.Z.) 
of the honeycomb lattice, indicating the position of the Dirac 
points. On the right, a three dimensional view of the energy 
spectrum shows the linear intersection of the bands at the 
two Dirac points. We denote the full bandwidth VK, and the 
minimum energy gap at the edges of the Brillouin zone Eg- 



quantum gases, a honeycomb lattice has recently been re- 
alised and investigated with a Bose-Einstein condensate 
[l6l[T7], but no signatures of Dirac points were observed. 
Here we study an ultracold Fermi gas of ^^K atoms in 
a two-dimensional tunable optical lattice, which can be 
continuously adjusted to create square, triangular, dimer 
and honeycomb structures. In the honeycomb lattice, we 
identify the presence of Dirac points in the band structure 
by observing a minimum band gap inside the Brillouin 
zone via interband transitions. Our method is closely re- 
lated to a technique recently used with bosonic atoms to 
characterize the linear crossing of two high-energy bands 
in a one-dimensional bichromatic lattice [18 , but pro- 
vides in addition momentum resolution. 

To create and manipulate Dirac points, we have devel- 
oped a two-dimensional optical lattice of adjustable ge- 
ometry. It is formed by three retro-reflected laser beams 
of wavelength A = 1064 nm, arranged as depicted in Fig. 
la. The interference of two perpendicular beams X and 
Y gives rise to a chequerboard lattice of spacing \/y/2. 
A third beam X, collinear with X but detuned by a fre- 
quency ^, creates an additional standing wave of spacing 
A/2. This yields a potential of the form 

V{x,y) = -V-^cos^{kx^e/2) - Vx^o^^ikx) 
—Vy cos^ {ky) — 2ol\JVxVy cos{kx) cos{ky) cos (p{l) 

where V^, Vx and Vy denote the single beam lattice 
depths (proportional to the laser beam intensities), a is 
the visibility of the interference pattern and k = 27r/A. 
We can adjust the two phases continuously, and choose 
= TT and (f = (see Methods). Varying the relative in- 
tensities of the beams allows us to realise various lattice 
structures, as displayed in Fig. lb. In the following we 
focus on the honeycomb lattice, whose real space poten- 
tial is shown in Fig. Ic. 

The honeycomb lattice consists of two sublattices A 
and B. Therefore, the wavefunctions are two-component 
spinors. Tunneling between the sublattices leads to the 
formation of two energy bands, which are well separated 
from the higher bands and have a conical intersection at 
two quasi-momentum points in the Brillouin zone - the 
Dirac points. These points are topological defects in the 
band structure, with an associated Berry phase of ±7r. 
This warrants their stability with respect to lattice per- 
turbations, so that a large range of lattice anisotropics 
only changes their position inside the Brillouin zone. In 
contrast, breaking the inversion symmetry of the poten- 
tial by introducing an energy offset A between sublattices 
opens an energy gap at the Dirac points, proportional to 
the offset. In our implementation, the sublattice offset 
A depends only on the value of the phase and can be 
precisely adjusted (see Methods). As displayed in Fig. 
Ic and d, the primitive lattice vectors are perpendicular, 
leading to a square Brillouin zone with two Dirac points 
inside. Their position is symmetric around the center and 
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FIG. 2: Probing the Dirac points, a, Quasi-momentum 
distribution of the atoms before and after one Bloch oscilla- 
tion of period Tb- The cloud explores several trajectories in 
quasi-momentum space simultaneously. For trajectory 1 (blue 
solid circle) the atoms remain in the first energy band. In con- 
trast, trajectory 2 (green open circle) passes through a Dirac 
point at t = Tb/2. There, the energy splitting between the 
bands vanishes and the atoms are transferred to the second 
band. When measuring the quasi-momentum distribution at 
t = Tb, these atoms are missing in the first Brillouin zone 
and appear in the second one. b, Dependence of the total 
fraction of atoms transferred to the second band ^ on the de- 
tuning of the lattice beams 6, which controls the sublattice 
energy offset A. The maximum indicates the point of inver- 
sion symmetry, where A = (^ = tt in eq. (1)) and the gap 
at the Dirac point vanishes. Away from the peak, the atoms 
behave as Dirac fermions with a tunable mass (see insets). 
Values and error bars denote the mean and standard devia- 
tion of five consecutive measurements, whereas the solid line 
is a Gaussian fit to the data. 
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fixed to Qx = 0^ owing to tiie time- reversal and reflection 
symmetries of the system [19j. 

We characterise the Dirac points by probing the en- 
ergy splitting between the two lowest energy bands 
through interband transitions [20^. The starting point 
of the experiment is a non-interacting ultracold gas of 
50,000 fermionic "^^K atoms in the iF^mp) = 
|9/2, — 9/2) state. The cloud is prepared in the lowest 
energy band of a honeycomb lattice with V-^ ^ y/^R — 
[4.0(2), 0.28(1), 1.8(1)], which also causes a 'weak har- 
monic conflnement with trapping frequencies uj^.y^zl"^^ = 
[17.6(1), 31.8(5), 32.7(5)] Hz. Here Er = h^/2mX^, h 
denotes the Planck constant and m the mass of ^^K. 
By applying a weak magnetic fleld gradient, the atomic 
cloud is subjected to a constant force F along the x di- 
rection, equivalent to an electric fleld in solid-state sys- 
tems. The atoms are hence accelerated such that their 
quasi-momentum Qx increases linearly up to the edge of 
the Brillouin zone, where a Bragg reflection occurs. The 
cloud eventually returns to the center of the band, per- 
forming one full Bloch oscillation [21 . We then measure 
the quasi-momentum distribution of the atoms in the dif- 
ferent bands [22 (see Methods). 

Owing to the flnite momentum width of the cloud, tra- 
jectories with different quasi- moment a Qy are simultane- 
ously explored during the Bloch cycle, as illustrated in 
Fig. 2a. For a trajectory far from the Dirac points, the 
atoms remain in the lowest energy band (trajectory 1). 
In contrast, when passing through a Dirac point (trajec- 
tory 2), the atoms are transferred from the flrst to the 
second band because of the vanishing energy splitting 
at the linear band crossing. When measuring the quasi- 
momentum distribution, these atoms are missing in the 
flrst Brillouin zone and appear in the second band, as can 
be seen in Fig. 2a. We identify the points of maximum 
transfer with the Dirac points. The energy resolution of 
the method is set by the characteristic energy of the ap- 
plied force [21^ Es/h = F\/2h = 88.6(7) Hz, which is 
small compared to the full bandwidth W/h = 4.6 kHz 
and the minimum band gap at the edges of the Brillouin 
zone Ec/h = 475 Hz. 

To investigate how breaking the inversion symmetry of 
the lattice affects the Dirac points, we vary the sublat- 
tice offset A, which is controlled by the frequency de- 
tuning S between the lattice beams, and measure the 
total fraction of atoms transferred to the second band 
^. The results obtained for a honeycomb lattice with 
^x,x,y/^R = [3.6(2), 0.28(1), 1.8(1)] are displayed in 
Fig. 2b, and show a sharp maximum in the transferred 
fraction. We identify this situation as the point of in- 
version symmetry A = (^ = tt), in good agreement 
with an independent calibration (see Methods). At this 
setting the band gap at the Dirac points vanishes. The 
population in the second band decreases symmetrically 
on both sides of the peak as the gap opens up, indicating 
the transition from massless to massive Dirac fermions. 




FIG. 3: Movement of the Dirac points, a, Distance of the 
Dirac points to the corners of the Brillouin zone, as measured 
through momentum-resolved interband transitions. The tun- 
neling along the x direction increases when decreasing the 
lattice depth V^. The distance is extracted from the second 
band quasi-momentum distribution after one Bloch cycle (see 
insets). The merging of the two Dirac points at the corners 
of the Brillouin zone is signalled by a single line of missing 
atoms in the first band. Values and error bars are the mean 
and standard deviation of three to nine measurements. The 
solid line is the prediction of a two-dimensional band structure 
calculation without any fitting parameters, b, Energy split- 
ting between the two lowest bands. It shows the displacement 
of the Dirac cones inside the Brillouin zone, as well as their 
deformation depending on the lattice depth Vj^. 



The relative strength of the tunnel couplings between 
the different sites of the lattice fixes the position of the 
Dirac points inside the Brillouin zone, as well as the slope 
of the associated linear dispersion relation [5-9 . How- 
ever, the tunability of our optical lattice structure allows 
for independent adjustment of the tunneling parameters 
along the x and y directions simply by controlling the 
intensity of the laser beams. For isotropic tunnelings the 
slope of the Dirac cones is the same in all directions, 
while being anisotropic otherwise. The distance of the 
Dirac points to the corners of the Brillouin zone along Qy 
can be varied between and (7b /2, whilst Qx = due to 
reflection symmetry [19]. Here qb = 27r/A denotes the 
Bloch wave vector. 

We exploit the momentum resolution of the inter- 
band transitions to directly observe the movement of the 
Dirac points. Starting from a honeycomb lattice with 
^x,x,y/^R = [5.4(3), 0.28(1), 1.8(1)], we gradually in- 
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crease the tunneling along the x direction by decreasing 
the intensity of the X beam. As displayed in Fig. 3, 
the position of the Dirac points continuously approaches 
the corners of the Brillouin zone, as expected from an 
ab initio two-dimensional band structure calculation (see 
Methods). 

When reaching the corners, the two Dirac points 
merge, annihilating each other. There, the dispersion 
relation becomes quadratic along the Qy axis, remain- 
ing linear along q^. Beyond this critical point, a finite 
band gap appears for all quasi-momenta of the Brillouin 
zone. This situation signals the transition between band 
structures of two different topologies, one containing two 
Dirac points and the other none. For two-dimensional 
honeycomb lattices at half filling, it corresponds to a Lif- 
shitz phase transition from a semi-metallic to a band- 
insulating phase [3 . 

We experimentally map out the topological transition 
line by recording the fraction of atoms transferred to the 
second band ^ as a function of the lattice depths 
and Vx, while keeping Vy/Er = 1.8(1). The results are 
displayed in Fig. 4a. There, the onset of population 
transfer to the second band signals the appearance of 
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FIG. 4: Topological transition. Fraction of atoms trans- 
ferred to the second band ^ as a function of lattice depths 
and Vx, with Vy/Er = 1.8(1). Different lattice ge- 
ometries (square, chequerboard, triangular, dimer and hon- 
eycomb) are realised, see Fig. lb. We consider trajectories 
in quasi-momentum space along the Qx and qy directions (a 
and b resp.). To maximise the transfer for the qy trajecto- 
ries, where the cloud successively passes the two Dirac points, 
we set = 1.013(l)7r. In both cases, the onset of population 
transfer to the second band signals the topological transition, 
where the Dirac points appear. The dashed line is the the- 
oretical prediction for the transition line without any fitting 
parameters. The insets show cuts of the band structure along 
the qy axis {qx = 0), illustrating the emergence of the Dirac 
points when increasing V^^. 



Dirac points in the band stucture of the lattice. The 
transferred fraction ^ decreases for large values of V^, as 
the Dirac point positions exceed the momentum width of 
the cloud. 

To extend the range of our measurements and probe 
the Dirac points even in this regime, we apply a force 
along the y direction. We hence explore a new class of 
trajectories in quasi-momentum space. This allows for 
the investigation of very anisotropic Dirac cones, which 
become almost ffat along the qx direction as we ap- 
proach the crossover to a one-dimensional lattice struc- 
ture {V^ ^ Vx)' Along the qy trajectories the center of 
the cloud successively passes the two Dirac points during 
the Bloch cycle, effectively realising a Stiickelberg inter- 
ferometer [23l [24] in a two-dimensional band structure. 
As displayed in Fig. 4b, we again identify the topological 
transition through the onset of population transfer to the 
second band. The results for the transition line obtained 
for both measurement series are in excellent agreement 
with ab initio band structure calculations. 

In this work we have realised Dirac points with highly 
tunable properties using ultracold fermionic atoms in a 
honeycomb optical lattice. A new class of physical phe- 
nomena is now within the domain of quantum gas exper- 
iments. Furthermore, the flexibility and novel observ- 
ables of these systems will provide new insights. For 
example, the unique coherence of quantum gases offers 
the possibility of directly measuring the Berry phase as- 
sociated to the Dirac points by interferometric meth- 
ods. Topological order could be achieved by introduc- 
ing artiflcial gauge flelds, either via Raman transitions 
[25 or time-dependent lattice modulation [26 . More- 
over, the exceptionally tunable lattice potential we have 
developed opens up a wealth of new avenues for opti- 
cal lattice experiments. For spin mixtures with repulsive 
interactions, the dynamic transition between dimer and 
square lattices should facilitate the adiabatic preparation 
of an anti-ferromagnetic phase [27 and enable the study 
of quantum criticality [28]. Additionally, the triangular 
and honeycomb lattices provide the possibility to explore 
magnetic frustration and spin liquid phases [Til US] • 

METHODS 

Preparation scheme. After sympathetic cooling 
with ^^Rb in a magnetic trap, 2 x 10^ fermionic ^^K 
atoms are transferred into a dipole trap operating at a 
wavelength of 826 nm. A balanced spin mixture of the 
|F,mF) = |9/2, -9/2) and |9/2,-7/2) states is evapo- 
ratively cooled at a magnetic fleld of 197.6(1) G, in the 
vicinity of the Feshbach resonance at 202.1 G. We obtain 
typical temperatures of 0.2 Tp, where Tp is the Fermi 
temperature. The fleld is subsequently reduced and a 
magnetic fleld gradient is used to selectively remove the 
|9/2,-7/2) component, while levitating the |9/2, -9/2) 
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atoms against gravity. This polarised Fermi gas is loaded 
into the two-dimensional optical lattice within 200 ms, 
and the dipole trap is then switched off. 

Tunable optical lattice. The optical lattice is pro- 
duced by the combination of three retro-reflected laser 
beams of linear polarization and wavelength A = 1064 
nm, red detuned with respect to the Di and D2 lines of 
^^K. The X and X beams propagate along the same axis, 
while the Y beam is at 90.0(1)°. The interference of the 
X and Y beams produces a chequerboard potential of 
the form [29^ Vi{x,y) = —Vxcos'^{kx) — VVcos^(%) — 
2a^/Vx\^ cos{kx) cos{ky) cos (p. The phase between the 
two beams at the position of the atoms is stabilised in- 
terferometrically to cp/Tv = 0.00(3) using a pair of ad- 
ditional beams detuned from each other and from X,X 
and Y. This results in a weak additional lattice along 
each axis of about O.lEji. The X beam creates a po- 
tential V2{x) = —V^cos^{kx + ^/2), where the phase 6 
determines the relative position of the chequerboard pat- 
tern and the one-dimensional standing wave. We control 
the value of in the center of the cloud by adjusting 
the frequency detuning 6 between the X and X beams. 
We infer the precise value 0/6 = [7r/384.7(6)] MHz"^ 
from the peak position in Fig. 2b. This is in good 
agreement with an independent calibration obtained us- 
ing Raman-Nath diffraction on a ^^Rb Bose-Einstein con- 
densate, which yields 0/S = [7r/388(4)] MHz-^ At 
the edges of the cloud the phase differs by approxi- 
mately ±10~^7r. The total lattice potential is given by 
Vi{x^y) + V2{x) and, depending on the relative inten- 
sities of the beams, gives rise to square, chequerboard, 
triangular and honeycomb lattices, as well as a staggered 
arrangement of dimers |30] and an array of weakly cou- 
pled one-dimensional chains. The visibility a = 0.90(5) 
and the lattice depths x y calibrated by Raman- 
Nath diffraction. The method has a systematic uncer- 
tainty of 10% for the lattice depths, whereas the statis- 
tical uncertainties are given in the main text. The two- 
dimensional lattice lies in the xy plane, whereas in the 
z direction the atoms are harmonically trapped. Owing 
to the absence of interactions the z direction decouples. 
The underlying trap frequencies in our system scale with 
the lattice depths according to the approximate expres- 
sions ujx (X >/TV, ^/V-^-\- {VxVy /Vx') and ujz oc 
Y^V^+ {VxVy /V^) + 1.241V. For the measurements in 
Fig. 2a we find a;^,^,, J27r = [17.6(1), 31.8(5), 32.7(5)] Hz, 
as calibrated from dipole oscillations of the cloud. 

Detection. The quasi-momentum distribution of the 
gas is probed using a band-mapping technique. The op- 
tical lattice beams are linearly ramped down in 500 /is, 
slowly enough for the atoms to stay adiabatically in their 
band while quasi-momentum is approximately conserved 
[22 . We then allow for 15 ms of ballistic expansion before 
taking an absorption image of the cloud. 

Band structure calculations. The energy spec- 
trum is obtained using an ab initio single-particle two- 



dimensional numerical band structure calculation for the 
homogeneous system. It therefore also takes into account 
higher order tunneling terms, which are relevant for the 
regime explored in this paper. In particular, they cause 
an asymmetry between the two lowest bands and lead to 
a tilt of the Dirac cones for certain parameter regimes. 
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